Nematic phase in the J1_1-J2_2 square lattice Ising model in an external field

The J$_1$-J$_2$ Ising model in the square lattice in the presence of an external field is studied by two approaches: the Cluster Variation Method (CVM) and Monte Carlo simulations. The use of the CVM in the square approximation leads to the presence of a new equilibrium phase, not previously reported for this model: an Ising-nematic phase, which shows orientational order but not positional order, between the known stripes and disordered phases. Suitable order parameters are defined and the phase diagram of the model is obtained. Monte Carlo simulations are in qualitative agreement with the CVM results, giving support to the presence of the new Ising-nematic phase. Phase diagrams in the temperature-external field plane are obtained for selected values of the parameter $\kappa=J_2/|J_1|$ which measures the relative strength of the competing interactions. From the CVM in the square approximation we obtain a line of second order transitions between the disordered and nematic phases, while the nematic-stripes phase transitions are found to be of first order. The Monte Carlo results suggest a line of second order nematic-disordered phase transitions in agreement with the CVM results. Regarding the stripes-nematic transitions, the present Monte Carlo results are not precise enough to reach definite conclusions about the nature of the transitions.Comment: 13 pages, 10 figure

Nematic phase in the J1-J2 square-lattice Ising model in an external field

© 2015 American Physical Society. The J1-J2 Ising model in the square lattice in the presence of an external field is studied by two approaches: the cluster variation method (CVM) and Monte Carlo simulations. The use of the CVM in the square approximation leads to the presence of a new equilibrium phase, not previously reported for this model: an Ising-nematic phase, which shows orientational order but not positional order, between the known stripes and disordered phases. Suitable order parameters are defined, and the phase diagram of the model is obtained. Monte Carlo simulations are in qualitative agreement with the CVM results, giving support to the presence of the new Ising-nematic phase. Phase diagrams in the temperature-external field plane are obtained for selected values of the parameter κ=J2/|J1| which measures the relative strength of the competing interactions. From the CVM in the square approximation we obtain a line of second order transitions between the disordered and nematic phases, while the nematic-stripes phase transitions are found to be of first order. The Monte Carlo results suggest a line of second order nematic-disordered phase transitions in agreement with the CVM results. Regarding the stripes-nematic transitions, the present Monte Carlo results are not precise enough to reach definite conclusions about the nature of the transitions.Peer Reviewe

Effects of confinement on pattern formation in two dimensional systems with competing interactions

Template-assisted pattern formation in monolayers of particles with competing short-range attraction and long-range repulsion interactions (SALR) is studied by Monte Carlo simulations in a simple generic model [N. G. Almarza et al., J. Chem. Phys., 2014, 140, 164708]. We focus on densities corresponding to formation of parallel stripes of particles and on monolayers laterally confined between straight parallel walls. We analyze both the morphology of the developed structures and the thermodynamic functions for broad ranges of temperature T and the separation L between the walls. At low temperature stripes parallel to the boundaries appear, with some corrugation when the distance between the walls does not match the bulk periodicity of the striped structure. The stripes integrity, however, is rarely broken for any L. This structural order is lost at T = T(L) depending on L according to a Kelvin-like equation. Above the Kelvin temperature T(L) many topological defects such as breaking or branching of the stripes appear, but a certain anisotropy in the orientation of the stripes persists. Finally, at high temperature and away from the walls, the system behaves as an isotropic fluid of elongated clusters of various lengths and with various numbers of branches. For L optimal for the stripe pattern the heat capacity as a function of temperature takes the maximum at T = T(L).Peer Reviewe

Phase behavior of the hard-sphere Maier-Saupe fluid under spatial confinement

The Maier-Saupe hard-sphere fluid is one of the simplest models that accounts for the isotropic-nematic transition characteristic of liquid crystal phases. At low temperatures the model is known to present a gas-liquid-like transition with a large difference between the densities of the coexistence phases, whereas at higher temperature the transition becomes a weak first-order transition resembling the typical order-disorder (nematic-isotropic) phase change of liquid crystals. Spatial dimensionality directly conditions the character of the orientational phase change (i.e., the high temperature transition), that goes from a first-order transition in the purely three-dimensional case, to a Berezinskii-Kosterlitz-Thouless-like continuous transition which occurs when the three dimensional Maier-Saupe spins are constrained to lie on a plane. In the latter instance, the ordered phase is not endowed with true long-range order. In this work we investigate how the continuous transition transforms into a true first-order phase change, by analyzing the phase behavior of a system of three dimensional Maier-Saupe hard spheres confined between two parallel plates, with separations ranging from the quasi-two-dimensional regime to the bulk three-dimensional limit. Our results indicate that spatial confinement in one direction induces the change from first order to a continuous transition with a corresponding decrease of the transition temperatures. As to the gas-liquid transition, the estimated “critical” temperatures and densities also decrease as the fluid is confined, in agreement with previous results for other simple systems.Dirección General de Investigación Científica y Técnica: MAT2007-65711-C04-04 Dirección General de Universidades e Investigación de la Comunidad de Madrid under Grant No. S0505/ESP/0299 and Program MOSSNOHO-CMPeer reviewe

Surface tension of the Widom-Rowlinson model

11 pags., 5 figs., 3 tabs.We consider the computation of the surface tension of the fluid-fluid interface for the Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] binary mixture from direct simulation of the inhomogeneous system. We make use of the standard mechanical route, in which the surface tension follows from the computation of the normal and tangential components of the pressure tensor of the system. In addition to the usual approach, which involves simulations of the inhomogeneous system in the canonical ensemble, we also consider the computation of the surface tension in an ensemble where the pressure perpendicular (normal) to the planar interface is kept fixed. Both approaches are seen to provide consistent values of the interfacial tension. The issue of the system-size dependence of the surface tension is addressed. In addition, simulations of the fluid-fluid coexistence properties of the mixture are performed in the semigrand canonical ensemble. Our results are compared with existing data of the Widom-Rowlinson mixture and are also examined in the light of the vapor-liquid equilibrium of the thermodynamically equivalent one-component penetrable sphere model. © 2007 American Institute of Physics.Financial support is due to the Spanish Dirección General de Investigación Project Nos. FIS2004-06627-C02-01 E.d.M. and FIS2004-02954-C03-01 N.G.A. and from Universidad de Huelva and Junta de Andalucía. Additional funding from the Dirección General de Universidades e Investigación Comunidad de Madrid, Spain under the MOSSNOHO-CM program Grant No. S0505/ESP/0299 and from the Engineering and Physical Sciences EPSRC of the UK Grant Nos. GR/N20317, GR/N03358, GR/N35991, GR/R09497, and EP/E016340, the Joint Research Equipment Initiative JREI GR/M94427, and the Royal Society Wolfson Foundation refurbishment grant is also acknowledged. Finally we are grateful to the Royal Society for the award of a International Short Visit grant which has facilitated the collaborative work

Demixing and confinement of non-additive hard-sphere mixtures in slit pores

© 2015 AIP Publishing LLC. Using Monte Carlo simulation, we study the influence of geometric confinement on demixing for a series of symmetric non-additive hard spheres mixtures confined in slit pores. We consider both a wide range of positive non-additivities and a series of pore widths, ranging from the pure two dimensional limit to a large pore width where results are close to the bulk three dimensional case. Critical parameters are extracted by means of finite size analysis. As a general trend, we find that for this particular case in which demixing is induced by volume effects, the critical demixing densities (and pressures) increase due to confinement between neutral walls, following the expected behavior for phase equilibria of systems confined by pure repulsive walls: i.e., confinement generally enhances miscibility. However, a non-monotonous dependence of the critical pressure and density with pore size is found for small non-additivities. In this latter case, it turns out that an otherwise stable bulk mixture can be unexpectedly forced to demix by simple geometric confinement when the pore width decreases down to approximately one and a half molecular diameters.Peer Reviewe

Theory of repulsive charged colloids in slit-pores

9 pags., 8 figs., 1 tab.Using classical density functional theory (DFT) we analyze the structure of the density profiles and solvation pressures of negatively charged colloids confined in slit pores. The considered model, which was already successfully employed to study a real colloidal (silica) suspension [S. H. L. Klapp, Phys. Rev. Lett. 100, 118303 (2008)10.1103/PhysRevLett.100.118303], involves only the macroions which interact via the effective Derjaguin-Landau-Verwey-Overbeek (DLVO) potential supplemented by a hard core interaction. The solvent enters implicitly via the screening length of the DLVO interaction. The free energy functional describing the colloidal suspension consists of a hard sphere contribution obtained from fundamental measure theory and a long range contribution which is treated using two types of approximations. One of them is the mean field approximation (MFA) and the remaining is based on Rosenfelds perturbative method for constructing the Helmholtz energy functional. These theoretical calculations are carried out at different bulk densities and wall separations to compare finally to grand canonical Monte Carlo simulations. We also consider the impact of charged walls. Our results show that the perturbative DFT method yields generally qualitatively consistent and, for some systems, also quantitatively reliable results. In MFA, on the other hand, the neglect of charge-induced correlations leads to a breakdown of this approach in a broad range of densities. © 2012 American Institute of Physics.A.G. and N.G.A. acknowledge financial support from the Dirección General de Investigación Científica y Técnica under Grant No. FIS2010-15502, and from the Dirección General de Universidades e Investigación de la Comunidad de Madrid under Grant No. S2009/ESP-1691 and Program MODELICO-CM. S.G. and S.H.L.K. would like to thank the German Research Foundation (DFG) for financial support within the Collaborative Research Center (CRC) 951 “Hybrid Inorganic/organic systems for Optoelectronics” and the International Graduate School (IRTG) 1524 “Selfassembled soft-matter nanostructures at interfaces.

Pattern formation in binary fluid mixtures induced by short-range competing interactions

Molecular dynamics simulations and integral equation calculations of a simple equimolar mixture of diatomic molecules and monomers interacting via attractive and repulsive short-range potentials show the existence of pattern formation (microheterogeneity), mostly due to depletion forces away from the demixing region. Effective site-site potentials extracted from the pair correlation functions using an inverse Monte Carlo approach and an integral equation inversion procedure exhibit the features characteristic of a short-range attractive and long-range repulsive potential. When charges are incorporated into the model, this becomes a coarse grained representation of a room temperature ionic liquid, and as expected, intermediate range order becomes more pronounced and stable

Phase diagram of a two-dimensional lattice gas model of a ramp system

Using Monte Carlo Simulation and fundamental measure theory we study the phase diagram of a two-dimensional lattice gas model with a nearest neighbor hard core exclusion and a next-to-nearest neighbors finite repulsive interaction. The model presents two competing ranges of interaction and, in common with many experimental systems, exhibits a low density solid phase, which melts back to the fluid phase upon compression. The theoretical approach is found to provide a qualitatively correct picture of the phase diagram of our model system.Comment: 14 pages, 8 figures, uses RevTex